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We investigate how the extragalactic proton component derived within the “escape model” can 
be explained by astrophysical sources. We consider as possible cosmic ray (CR) sources nor- 
mal/starburst galaxies and radio-loud active galactic nuclei (AGN). We find that the contribution 
to the total extragalactic proton flux from normal and starburst galaxies is only subdominant and 
does not fit the spectral shape deduced in the escape model. In the case of radio-loud AGN, we show 
that the complete extragalactic proton spectrum can be explained by a single source population, 

BL Lac/FR I, for any of the potential acceleration sites in these sources. We calculate the diffuse 
neutrino and 7 -ray fluxes produced by these CR protons interacting with gas inside their sources. 

For a spectral slope of CRs close to a = 2.1 — 2.2 as suggested by shock acceleration, we find that 
these UHECR sources contribute the dominant fraction of both the isotropic 7 -ray background and 
of the extragalactic part of the astrophysical neutrino signal observed by IceCube. 


I. INTRODUCTION 

The search for the sources of ultrahigh energy cosmic 
rays (UHECR) and for an understanding of their acceler¬ 
ation mechanism is one of the important challenges in as- 
troparticle physics. It has been hoped for that CRs at the 
highest energies would be only weakly deflected in mag¬ 
netic fields and hence sources could be identified by usual 
astronomical methods. Both the Pierre Auger Observa¬ 
tory (PAO) and the Telescope Array (TA) do observe 
anisotropies in the arrival directions of UHECRs 01 
and a hot spot observed by TA has also a large statisti¬ 
cal significance. However, these anisotropies extend over 
medium angular scales, similar as found in Ref. Q for 
previous experimental data, and no successful correla¬ 
tion of UHECRs with potential astrophysical sources has 
been achieved yet. Meanwhile, there has been steady 
progress in other areas: The all-particle CR spectrum 
has been measured precisely, and data on the primary 
composition have become available. The Auger collabo¬ 
ration provided in 1 constraints on the fraction of four 
different elemental groups above 6 x 10 17 eV, while the 
KASCADE-Grande experiment covered with its compo¬ 
sition measurements energies up to 2 x 10 17 eV 0. Al¬ 
though the derived composition depends on the hadronic 
interaction models used for the analysis, the following 
qualitative conclusions can be drawn: First, the proton 
fraction amounts to ~ 40-60% in the energy range be¬ 
tween 7 x 10 17 eV and 7 x 10 18 eV and decreases after¬ 
wards, while the fraction of intermediate nuclei (He, N) 
increases. Second, the iron fraction in the energy range 
between 7 x 10 17 eV and 2 x 10 19 eV is limited by < 15- 
20% and its central value is consistent with zero. Thus 
the Galactic contribution to the observed CR spectrum 
has to die out around 7 x 10 17 eV, unless an additional 
sub-dominant and heavy Galactic component remains. 


Both light and intermediate elements above this energy 
cannot have a Galactic origin because their anisotropy 
would otherwise overshoot the upper limits set by Auger 
on the CR dipole anisotropy, see Ref. Q. 

These results provide a strong constraint on models for 
the transition between the Galactic and the extragalac¬ 
tic CR component. In particular, they exclude the dip 
model (7] which requires a proton fraction > 90%. In 
other models, the ankle at E ss a few x 10 18 eV has been 
identified with the transition between Galactic and ex¬ 
tragalactic CRs Si- However, such a high value of the 
transition energy exaggerates the acceleration problem 
of Galactic CR sources and contradicts the low iron frac¬ 
tion determined by Auger. Alternatively, it has been sug¬ 
gested that two populations of extragalactic CR sources 
exist, one dominating below and one above the ankle E3- 
Since there exists no convincing model for these two 
source classes and their properties, a more economical 
explanation based on a single extragalactic source type 
is desirable. 

Several recent studies tried to explain the measured 
spectrum and the composition both above and below the 
ankle using models with a mixed composition 00. 
Two of them, Refs. 0 m, explain the light composi¬ 
tion below the ankle by photo-disintegration in the CR 
sources, while at higher energies the initial mixed com¬ 
position survives. 

In Ref. m, the energy spectra of nuclei in the source 
are assumed to follow a 1 /E power-law. As a result, the 
extragalactic sources do not contribute at E < 3 x 10 17 
eV. In Ref. 0, the original index of nuclei spectra is 
around 2.1 — 2.2, as required by acceleration models, 
but all low energy nuclei are photo-disintegrated in the 
sources. As a result, the spectra of nuclei leaving the 
source are close to a 1/E power-law, while the proton 
spectrum follows the original acceleration spectrum be- 
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cause of a decaying neutron component which escapes 
from the source. These models predict a contribution 
of extragalactic sources to the proton spectrum in the 
energy range measured by KASCADE-Grande. 

A disadvantage of this type of models is that they 
predict a negligible contribution from extragalactic CR 
sources to the observed IceCube neutrino signal. Thus 
the extragalactic neutrino flux is disconnected from the 
extragalactic CR flux, requiring that extragalactic neu¬ 
trinos are produced in “hidden sources” with a large in¬ 
teraction depth for protons. In both models, sources with 
strongly positive evolution give the major contibution to 
the UHECR flux. In contrast, the authors of Ref. 0 
study sources with negative evolution in order to explain 
the composition measured by Auger with softer fluxes 
oc 1/E 2 without invoking photo-disintegration. 


In this work, we investigate which CR source classes 
can explain the extragalactic proton component derived 
within the escape model 0 , 0 ] and, at the same time, 
can give a significant contribution to the isotropic 7 - 
ray background (IGRB) measured by Fermi-LAT [16| 
and to the astrophysical neutrino signal observed by 
IceCube 0 - These neutrinos have been suggested to 
have either a Galactic (e.g. 0-H3) and/or an extra¬ 
galactic origin (e.g. 0), and may have a hadronuclear 


origin [22M24I ]. Amongst others, radio-loud AGNs and 
hosting galaxy clusters have being proposed as possible 
sources 


221, and a connection to 10 17 ' 5-185 eV CRs was 
suggested and investigated in 0. Secondary neutrinos 
and 7 -rays from UHECR propagation in clusters were 
calculated in Ref. 1261. 


We consider here as possible CR sources nor¬ 
mal/star burst galaxies and radio-loud active galactic nu¬ 
clei (AGN). This choice is motivated by the fact that 
these sources might give the dominant contribution to 
the IGRB at low and high energies, respectively. In par¬ 
ticular, it has been suggested in Refs. 0] that blazars 
can contribute up to 100% to the IGRB. As production 
mechanism of the secondary 7 -ray and neutrino fluxes, 
we use CR interactions with gas in their sources. We 
find that normal and star-forming galaxies can explain 
neither the spectral shape nor the magnitude of the de¬ 
rived extragalactic proton flux. In the case of radio-loud 
AGNs, we show that the complete extragalactic proton 
spectrum can be explained by a single source popula¬ 
tion. We calculate also the diffuse neutrino and 7 -ray 
fluxes produced by these CR protons interacting with 
gas. For a spectral slope of CRs close to a p = 2.1 — 2.2 
as suggested by shock acceleration, we find that these 
UHECR sources can contribute the dominant fraction to 
the isotropic 7 -ray background (IGRB), and the major 
contribution to the extragalactic part of the astrophysical 
neutrino signal observed by IceCube. 



FIG. 1: Extragalactic proton flux deduced in Ref. El within 
the escape model from KASCADE-Grande and Auger data 
(brown errorbars) together with the original KASCADE and 
KASCADE-Grande proton data and the total CR flux from 
Auger. Also shown is the predicted proton flux from BL Lacs 
(red line). 


II. METHODOLOGY 


A. Extragalactic CR proton flux 


We summarize first how the extragalactic CR pro¬ 
ton flux in the escape model was derived in Ref. Il4| . 
In a first step, we derived the Galactic all-particle CR 
flux summing up all CR groups obtained for the max¬ 
imal rigidity lZ max = 10 17 V and accounting for the 
Auger iron constraint. Then we obtained the total ex¬ 
tragalactic flux by subtracting the predicted total Galac¬ 
tic flux from the measured total CR flux. The extra¬ 
galactic proton flux observed by Auger followed apply¬ 
ing the composition measurement [ 3 ] where we chose the 
results obtained using the EPOS-LHC simulation. Fi¬ 
nally, we deduced the contribution of extragalactic pro¬ 
tons to the observed proton flux by KASCADE and 
KASCADE-Grande, subtracting from the proton fluxes 
given in Ref. I5j the prediction of the escape model: Since 
the predicted Galactic proton flux lies for energies above 
E > 3 x 10 16 eV below the measured one, the difference 
has to be accounted for by extragalactic protons. We 
show in Fig. Q] the deduced extragalactic proton fluxes 
from KASKADE-Grande and Auger data in the escape 
model with brown errorbars. Combining the KASCADE, 
KASCADE-Grande and Auger data suggests that the 
slope of the extragalactic proton energy spectrum is flat 
at low energies, E < 10 18 eV, consistent with a p ~ 2.2, 
and softens to a p ~ 3 at higher energies, E > 10 18 eV, 
cf. Fig. [ 6 ] 
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B. CR interactions with gas and photons 

We summarize now how we calculate the interactions 
of CRs with gas and the extragalactic background light 
(EBL). We split the propagation in two parts: The first 
one includes the propagation in the source, the host 
galaxy and galaxy cluster where we assume that pro¬ 
ton interactions with gas are dominant. The spectrum 
of exiting particles is then used in the second step as 
an “effective source spectrum” from which we calculate 
the resulting diffuse flux taking into account the distri¬ 
bution p(z, L ) of sources as well as the interaction of pro¬ 
tons, electron and photon with the EBL and the CMB. 
For both steps, we use the open source code [28} which 
solves the corresponding kinetic equations in one dimen¬ 
sion. We employ the baseline EBL model of Ref. [29[. 

As input for the first step, we require the energy depen¬ 
dent grammage X(E) and the proton injection spectrum 
dNcFt/dE. Starting from the injection spectrum of pro¬ 
tons, we simulate their propagation and their interaction 
to obtain the spectra of protons and secondary particles 
leaving the “effective source”. We neglect all interactions 
except pp interactions in the “effective source”. This as¬ 
sumption is not satisfied for some UHECR acceleration 
sites in the vicinity of the AGN, such as notably acceler¬ 
ation close to the accretion disk where the intense radia¬ 
tion field would make py interactions dominate over pp. 
However, for acceleration at the polar caps, the radiation 
field is sufficiently low for pp to dominate. We also as¬ 
sume that all secondaries escape freely, except electrons. 
The fate of electrons depends on the strength of the ra¬ 
diation and magnetic fields and on the source size. In ex¬ 
tended sources as galaxies with relatively small magnetic 
fields, they lose all their energy via synchrotron and in¬ 
verse Compton radiation. For simplicity, we neglect pair 
production by photons inside the source, since the fol¬ 
lowing cascade development outside the source leads to 
a universal spectrum. 

The code [H} has been extended implementing pp in¬ 
teractions as follows: The inelastic cross sections CTi ne i 
of CR nuclei on gas were calculated with QGSJET-II- 
04 [30| . For the spectrum of secondary photons and neu¬ 
trinos produced in pp interactions differential cross sec¬ 
tions tabulated from QGSJET-II-04 were used [Tl! . Sec¬ 
ondaries from heavier elements in the CR flux are sup¬ 
pressed. For this reason, we will only need to consider the 
extragalactic CR proton flux in the following, when try¬ 
ing to fit CR, gamma-ray and neutrino fluxes altogether. 
The contributions from CR nuclei are included adding 
a nuclear enhancement factor Em- Similarly, the helium 
component of the interstellar medium was accounted for. 
Combining both effects, we set Em = 2.0 [H|. In order to 
take into account properly the energy dependence of the 
grammage, while still using kinetic equations in a one¬ 
dimensional framework, we include in the pp interaction 
rates R{E) the energy dependent grammage X(E). In 
the case of radio-loud AGNs, we have used for the gram¬ 
mage the simple parameterization X(E) oc E _1 / 3 ex¬ 


pected for a turbulent magnetic field with a Kolmogorov 
spectrum. The normalization was fixed by setting the in¬ 
teraction depth t pp = 1 at a reference energy E esc , which 
is the only free model parameter in this case. In contrast, 
for star-forming galaxies, we have used the grammage de¬ 
rived within the ecape model, see next section. 


III. STAR-FORMING GALAXIES 

The escape model developed in Refs. 00 provides 
an excellent description of the measured fluxes of CR nu¬ 
clei from below the knee to ~ 10 18 eV. The model predicts 
an early transition from Galactic to extragalactic CRs. 
Natural candidate sources for (part of) the extragalactic 
flux are all other star-forming galaxies (i.e. normal spi¬ 
ral galaxies, starburst galaxies and star-forming AGNs). 
Moreover, starburst galaxies are expected to give a ma¬ 
jor contribution to the IGRB, and they have also been 
considered as possible sources of high-energy neutrinos, 
assuming they are CR calorimeters -see amongst oth¬ 
ers Refs. 1 3 (it . Therefore, we first examine how large 
the contribution from star-forming (SF) galaxies to these 
fluxes is. 


A. CR flux from a single galaxy 

We apply the escape model to other normal spiral 
galaxies. In Refs. [lj,[l5|, we argued that a Kolmogorov 
spectrum for the turbulent galactic magnetic field to¬ 
gether with a 1/E 22 power-law injection spectrum pro¬ 
vides an explanation of the knee and of all CR compo¬ 
sition measurements between 10 14 eV and 10 18 eV. The 
deviation of the local proton spectrum from this slope 
is naturally explained by the influence of a local, recent 
source [37}. Therefore we use a universal 1/E 2 - 2 power- 
law as injection spectrum for all nuclei. Inside galaxies, 
the spectrum is modified because of energy-dependent 
CR confinement in them, whereas the CR flux exiting 
them retains the original injection 1/E 2 - 2 spectrum. We 
assume that spiral galaxies have magnetic fields which are 
on average similar to the one in the Milky Way. Hence we 
can use for the confinement time t(E) oc X(E ) of CRs 
in norm al sp iral galaxies the grammage X (E) calculated 
in Refs, [lj, [l5| for the Milky Way. 

In contrast, the observed strengths of magnetic fields 
in starburst galaxies are a factor ~ 100 larger than in 
the Milky Way [38} . Since the confinement time t(E) 
is a function of E/(ZeB), we can compensate for this 
change in B by rescaling the CR energy E in the gram¬ 
mage calculated for the Milky Way. We assume that the 
coherence length l c of the turbulent interstellar magnetic 
fields inside starburst galaxies is not significantly differ¬ 
ent from that in normal galaxies. 

In Fig. [2] we show the resulting interaction depth 
r = ai ne iX/mp in a starburst galaxy for various CR nu¬ 
clei as a function of energy, where we used the density 
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FIG. 2: Interaction depth r in a starburst galaxy, as a func¬ 
tion of CR energy E, for Fe (top, green line), A1 (magenta), 
C (black), He (blue) and protons (bottom, red). 



z 


FIG. 3: Rate p of type II SNe (green curve) and effective 
comoving source density N c (z) of BL Lacs (red) as a function 
of redshift, normalized to their present values at z = 0. 


profile described in m with no = 1/cm 3 . The fraction of 
CRs which are absorbed is given by /i n t = 1 — exp(— t). 
Thus the flux of heavy nuclei, such as A1 or Fe, leav¬ 
ing a starburst galaxy is exponentially suppressed. In 
contrast, CR interactions on gas can be neglected in nor¬ 
mal galaxies because of their weaker magnetic fields and 
shorter CR confinement times. 

Let us assume that every normal spiral galaxy (resp. 
every starburst galaxy) accelerates CRs with charge Z 
up to a maximum energy E max = Z x 10 17 eV (resp. 
E max = Zxl0 18 eV). The larger maximum energy achiev¬ 
able in starburst galaxies may e.g. be connected to the 
larger magnetic fields, which makes an additional acceler¬ 
ation of CRs in superbubbles more efficient [39|. Finally, 
we assume that the fractions of the injected CR nuclei are 
the same as the ones below the knee in the Milky Way, 
except for the proton fraction. We assume that the local 
proton fraction measured just below the knee is reduced 
due to the steeper proton spectrum compared to the av¬ 
erage one in normal galaxies. We choose to set the proton 
fraction to f p = 0.5 of the CR flux. Thus the composi¬ 
tion of p:He:N:Al:Fe is fixed as 50:22.5:12.5:5:10, which 
is consistent with the one deduced in 0 for He:N:Al:Fe 
at E = 10 14 eV. 


B. Diffuse fluxes from normal spiral and starburst 
galaxies 

We compute now the diffuse flux of extragalactic CRs 
from all star-forming galaxies. Let us first assume that 
at redshift z the CR emissivity Qcr(z) scales with the 
Type II supernova (SN) rate psnii(^), and thus with the 
star-formation rate (SFR). This is a phenomenological 
argument, which does not directly imply that these high- 
energy CRs have been accelerated at supernova shock 
waves -see a for a discussion on possible acceleration 
mechanisms. For /OsNii(-t): we use the parametrization 


presented in [4(| and take the Baldry & Glazebrook ini¬ 
tial mass function |0]. The corresponding Type II SN 
rate is [40j 


Psnii(z) 


0.0132 (a + bz ) h 
1 + ( z/c) d 


yr _1 Mpc -3 , 


(1) 


with a = 0.0118, b = 0.08, c = 3.3, d = 5.2 and h = 0.7. 
It is plotted in Fig. [3] as a green line, which peaks at 
redshift z ~ 2 — 3. We assume that, globally, a fraction 
£CR = 0.1 of the kinetic energy of supernovae (-Esn ~ 
10 51 erg per SN) is channeled into CRs. The integral CR 
emissivity <3cr(~) is then 


Qgr(z) ~ 9 x IQ’ 22 


(a + bz) h 
1 + ( z/c) d 


eV 
cm 3 s 


( 2 ) 


We define Qcr.SB (resp. Qcr.Sp) as the integral CR 
emissivity due to starburst (resp. normal spiral) galaxies. 
Qcr,SB = /sb Qcr and Qcr,SP = (1 - /sb) Qcr, where 
/sb denotes the fraction of SFR (or SNe) occurring in 
starburst galaxies at redshift z. We parametrize /sb with 
the two examples proposed in 0 : (*) /sb = 0.9 z + 0.1 
at z < 1, and = 1 otherwise, or (ii) /sb = 0.1 (1 + z) 3 at 
z < 1, and = 0.8 otherwise. We also consider a scenario 
(in) motivated by Ref. [H|, where star-forming galax¬ 
ies are divided into four classes: normal spiral galaxies, 
starburst galaxies and star-forming AGNs. The latter 
category is divided into two subsets: SF-AGNs which 
resemble spiral galaxies and SF-AGNs which resemble 
starburst galaxies. 

We can now determine the differential CR emissivity 
9cr,Sb(E, z) from starburst galaxies, with 






(3) 

taking an injection spectrum ^ oc E 2 ' 2 , between 
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-Emin = 1 GeV and -E max = 10 18 eV. This yields 


qcR,i{E , z) 


1 x 10~ 2 ° (o + bz) h ( E \ 
eVcm 3 s 1 + {z/c) d \leVy 


with i = {SB, SP}. 

We can find the resulting diffuse CR intensity 1(E) 
from 


1(E) 


c 

4t tH 0 



dz 

(1 + z)u) 


9CR (z,E') 


e ~ T ( E ') , 


(5) 


where ui = + U m (l + z) 3 . In the following, we take 

H 0 = 70kms _1 Mpc -1 , Qa = 0.7 and O m = 0.3. Be¬ 
low 10 18 eV, we can neglect CR interactions in the inter- 
galactic space and thus set E' = (1 + z)E. For such CR 
energies, we can neglect absorption during propagation, 
and t corresponds to the one given by Fig. [2] The upper 
integration limit z max is given by the magnetic horizon. 
In the following, we consider the optimistic case where 
the intergalactic magnetic fields are sufficiently weak to 
have z max > 6 for E > 10 16 eV. This case of no magnetic 
horizon yields an upper limit on the diffuse extragalactic 
CR flux one can expect from all star-forming galaxies. 
We set psnii(z) to zero at redshifts z > 6, as in lioj . 
Then, 


Isb(E) 


c 

AttHq 


10 - 2 ° 
eV cm 3 s 



x 



— fsB (1 + Z) 

LU 


- 3.2 


(a + bz) h 
1 + (z/c) d ' 


One finds for the total diffuse CR intensity 1(E ), 


E 2 1(E) 


1 x 10 2 eV / E \~°' 2 
cm 2 s sr \ 10 15 eV J 


( 6 ) 


and I S b(E) ~ 0.63 1(E) (resp. I SB (E) ~ 0.46/(.E)) in 
scenario (*) (resp. scenario (m)). With f p = 0.5, this 
gives a CR proton flux that is more than one order of 
magnitude weaker than the extragalactic CR proton flux 
deduced for the escape model: See the orange line in 
Figure 11 of Ref. [l4[. This shows that the guaranteed 
contribution from all star-forming galaxies in the escape 
model to the diffuse CR flux is negligible. Note that this 
conclusion does not depend on our assumption about the 
maximal acceleration energy in star-forming and star- 
burst galaxies. 

Using the above parameters, we present in Fig. [I] (left 
panel) the diffuse fluxes of CR nuclei due to normal and 
starburst galaxies between E = 10 16 eV and E = 10 18 eV. 
In this computation, we use scenario (*) for the evolution 
with redshift of the fraction of SNe occurring in starburst 
galaxies. We consider IGMFs with strength B = 10 -17 G 
and coherence length l c = 1 Mpc, e.g. values which are 
consistent with lower limits from y-ray observations for 
time-varying sources [44| . With such parameters, there 


is no magnetic horizon for any of the nuclei at the ener¬ 
gies we consider. The individual contributions of protons 
and nuclei from normal spiral (resp. starburst) galaxies 
are shown with the dashed (resp. solid) lines below the 
total flux. Orange lines for protons, blue ones for he¬ 
lium, green ones for CNO nuclei, red ones for aluminium 
and magenta ones for iron. The fluxes of intermediate 
and heavy nuclei from starburst galaxies are exponen¬ 
tially suppressed at the lowest energies, because of nuclei 
suffering significant energy losses on background gas in 
starbursts, see Fig. [2] In Fig. [I] (right panel), we show, 
for comparison, the results for scenario (in), where the 
CR flux is dominated by SF-AGNs. 


We consider next the contribution from star-forming 
galaxies to the diffuse neutrino and y-ray fluxes. The 
fraction of CR protons that interact in starburst galaxies 
is given by /; nt = 1 — exp(—r), where r is the interac¬ 
tion depth for protons shown in Fig. [2] We present in 
Fig. [5] the resulting gamma-ray (red line) and neutrino 
(magenta) flux from starburst galaxies, within evolution 
scenario ( i ). We also plot the diffuse CR proton flux 
from starburst galaxies. First, we note that star-forming 
galaxies give a subdominant contribution to the primary 
CR flux, for the spectral index a p = 2.2 which is favored 
by the escape model. This discrepancy could be reduced 
by increasing e.g. the fraction £cr of energy transferred 
to CRs or the SN rate. However, the redshift evolu¬ 
tion of star-forming galaxies leads to a proton flux which 
spectral shape disagrees with the shape deduced for the 
extragalactic proton flux. Thus we conclude that star¬ 
forming galaxies cannot be the main contributor to the 
extragalactic proton flux (up to their E max ). Next, we 
compute the secondary neutrino and photon fluxes. We 
compare them respectively to the astrophysical neutrino 
flux (magenta errorbars) measured by IceCube E3 and 
to the measurement (red errorbars) of the extragalactic 
y-ray background (EGB) by Fermi-LAT [Ij|. For the lat¬ 
ter, we show both the IGRB (lower curve) and the total 
EGB including resolved sources (upper curve). Choosing 
as average gas density n = 1/cm 3 , star-forming galax¬ 
ies contribute around 30% at 10 GeV to the IGRB, while 
their contribution to the neutrino signal observed by Ice- 
Cube reaches 10% below 10 14 eV. Increasing the gas den¬ 
sity by a factor of ten leads only to an increase of the 
secondary fluxes by a factor a few, because the source 
is already thick in a large energy range. Taking into ac¬ 
count the uncertainty in the grammage used, we conclude 
that star-forming galaxies cannot explain the extragalac¬ 
tic proton flux, but may contribute a significant fraction 
to the IGRB and to the extragalactic part of the neutrino 
signal as observed by IceCube, especially at low energies. 


Not surprisingly, studies which assume a p ~ 2.0 and 
a large grammage in starburst galaxies can reproduce a 
high-energy neutrino flux comparable to that of IceCube, 
see e.g. jMlH- 
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FIG. 4: Left panel: Diffuse CR flux from normal spiral (dashed lines) and starburst galaxies (solid lines) within scenario (i), 
as a function of energy. Total flux in black. Individual contributions of protons (orange), He (blue), N (green), A1 (red) 
and Fe (magenta). Right panel: Diffuse CR flux from normal spirals (dashed lines), star-forming AGNs -both starburst and 
non-starburst ones- (thick solid line) and starburst galaxies (thin solid lines) as a function of energy, within scenario (iii). For 
both panels, IGMF with strength B = 10” 17 G and coherence length l c = 1 Mpc. 



FIG. 5: Diffuse CR proton (blue), gamma-ray (red) and neu¬ 
trino (magenta) fluxes from starburst galaxies together with 
CR proton data from KASCADE, KASCADE-Grande 0 
and Auger (black errorbars) 0 , ll-the IGRB and the EGB 
from Fermi-LAT (red errorbars) [lg], and hi gh-e nergy neu¬ 
trino flux from IceCube (magenta errorbars) [131. Solid red 
and magenta lines for gamma-ray and neutrino fluxes with 
n = lcm -3 , and dashed lines for n = 10cm” 3 . 


IV. UHECR SOURCES 


the one of star-forming galaxies. Second, BL Lacs have 
been su gges ted to be a major contributor to the IGRB 
(see e.g. [46H48I ]). which raises the question whether these 
sources can also fit at the same time the extragalactic CR 
proton flux expected in the escape model. Also, AGNs 
are natural candidates of high-energy neutrino sources, 
see amongst others Refs. [49l. [50|. and Ref. [51!] for neu¬ 
trino production in the inner jets of radio-loud AGNs, 
including blazars, as sources of UHECRs. 

A. Evolution of BL Lacs 

We determine the cosmological evolution of BL 
Lac/FR I sources from the corresponding evolution of 
the 7 -ray luminosity, assuming that the CR and the 7 - 
ray luminosity are proportional, 

jjmax 

N c (z) oc / p(z, L 7 )L 1 cLL 1 . (7) 

J L,min 

7 

Here, p(z,L~/) is the 7 -ray luminosity function (LF), i.e. 
the number of sources per comoving volume and lumi¬ 
nosity. For p(z,Lry) we adopt the Luminosity-Dependent 
Density Evolution (LDDE) model of Ref. [47]. Within 
this model, the LF p(z,L 7 ) can be expressed as 


We discuss next AGNs, which are generally considered 
to be prime candidates for the sources of UHECRs. We 
consider the subset of radio-loud AGN or, more precisely, 
the BL Lac/Fanaroff-Riley I (FRI) sub-class of the radio- 
loud AGN population. This choice is motivated by two 
reasons: First, the evolution of these sources is relatively 
slow and peaks at low redshift. Thus, the resulting dif¬ 
fuse CR flux has a rather different spectral shape than 


p(^,i 7 ) = p(L 7 )e( 2 :,i 7 ), ( 8 ) 


with 


p(Lj) — 

e ( Z ’ L-y) = 


A 


log ( 10 )L 7 
l + z 
_ V1 + z c{L 7 ) 






l + z 


1 + z c (L 7 ) 


, (9) 


( 10 ) 
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and 


Zc{ - L7) Z * (lO 48 ergs -1 ) ' (U) 

The numerical values for the parameters were determined 
in [47] from a fit to the statistics of BL Lacs observed by 
the Fermi-LAT telescope, cf. their Table 3. The evolution 
of the effective source density with the redshift is shown 
in Fig. [3] In contrast to average AGNs, the number den¬ 
sity of BL Lac and FR I galaxies peaks at low redshift, 
z < 1. Their evolution is similar to that of galaxy clus¬ 
ters. In fact, most of the FR I sources, which are in the 
unified AGN scheme the same sources as BL Lacs seen 
under different observation angles, reside in the centres of 
the dominant central elliptical galaxies of galaxy clusters 
(cD galaxies). 


velocity is also v ~ c because the jet is a relativistic out¬ 
flow. Only in the case of the large scale radio lobes the 
velocity could be r < c. In this case it is determined by 
the details of interaction of the lobes with the interstel¬ 
lar/intracluster medium. In any case, the eddy turnover 
scale is certainly much shorter than the source lifetime 
for all the elements of the radio loud AGN. This means 
that the medium and magnetic field in the source are 
turbulent. 

The turbulence power spectrum may for example fol¬ 
low a Kolmogorov or an Iroshnikov-Kraichnan power law. 
In our calculations, we assume that the power spectrum 
of the turbulence is a Kolmogorov one. Then the escape 
time of CRs scales with energy as 


R / E \ ^ 1/3 

C \ -Efree J 


5xl0 6 s 


E 

10 11 eV 



(14) 


B. Interactions in BL Lac/FR I sources 

We assume that the CR injection spectrum of each 
source follows a power-law with slope a p and exponential 
cut-off, 


dNpR 

dE 


oc E ap exp 



( 12 ) 


For each assumed slope a p of the spectrum, we adjust 
the cutoff energy E cut in such a way that the spectrum 
of the entire source population (integrated over redshift) 
fits best the observed cosmic ray spectrum in the energy 
range 10 17 eV 10 20 eV. 

Cosmic rays of low energy are not necessarily escap¬ 
ing from the source. First of all, they could be trapped 
right in the source. The condition of free escape from the 
source is that the Larmor radius of the accelerated par¬ 
ticle is comparable to the source size R. This condition 
reads E > Ef ree , where 


Kf ree ~ eBR ~ 3 x 10 2 ° eV 


B 

10 4 G 


R 

10 14 cm ’ 


(13) 


e is electric charge of the particle and B is the magnetic 
field strength. Lower energy particles are trapped inside 
the source, be it the AGN central engine, jet or the radio 
lobes. 

The trapped particles can still escape from the source, 
but in a diffusive way on a much longer time scale. The 
details of this process depend on the turbulence of the 
magnetic field in the relevant source structure. The time 
scale of turbulence development on a distance scale A 
can be estimated by the eddy turnover time T tur b ~ X/v, 
where v is the average bulk velocity of the plasma moving 
over the distances of the order of A. In the case of the 
central engine of AGN, this velocity scale is the typical 
velocity of the accretion flow, which is about the Keple- 
rian or free-fall velocity. Close to the black hole horizon, 
this velocity is relativistic, v ~ c. In the AGN jet, the 


where we used as source size R = 10 14 cm and as mag¬ 
netic field strength B = 10 4 G. 

Cosmic rays trapped inside the source lose energy 
by interacting with the ambient medium present in the 
source. In the case of the accretion flow, the density is 
moderately low n < 10 10 cm -3 for the radiatively ineffi¬ 
cient accretion flows powering FR I/BL Lac sources. The 
energy loss time of CR protons is 

tpp = cna pp n “ 1 x 10 s (i 0 9 cm" 3 ) ’ 

where a pp ~ (3 — 8) x 10 ~ 26 cm 2 is the inelastic pp cross- 
section and k — 0.6 the inelasticity. The interaction time 
is shorter than the escape time, t esc > t pp , for CRs with 
energy E < E esc = 8 x 10 12 eV, where we used again R = 
10 14 cm and B = 10 4 G for the numerical estimate. Thus 
CRs with energies below ~ 10 TeV would not escape from 
the central engine of an AGN powered by a 3 x 10 8 Mq 
black hole. Note however that the numerical value of the 
escape energy depends strongly on the chosen values for 
n, B and R and should be considered therefore only as 
an indication. 

Our results for the diffuse flux of CR protons from 
UHECR sources following the BL Lac evolution ([9j) is 
shown in Fig. [0] for E cut = 10 19 eV, and two different 
values of CR slope, a p = 2.17 for the upper panel and 
a p = 2.10 for the middle and lower panels. The Galactic 
proton flux in the escape model is shown with a dashed 
blue line. The choice a p = 2.17 (resp. a p = 2.1) results 
in an excellent (resp. good) fit of the extragalactic proton 
component deduced from Auger and KASCADE-Grande 
measurements. 

In the same figure, we also show the secondary fluxes 
obtained for E esc = 3 x 10 11 eV (upper and middle pan¬ 
els) and Eesc = 10 14 eV (lower panel), i.e. values of E esc 
which are characteristic for CR acceleration close to the 
supermassive BH powering the BL Lac. The diffuse pho¬ 
ton (resp. neutrino) fluxes are shown with red (resp. 
magenta) lines. 
















FIG. 6: Diffuse flux of CR protons from BL Lacs (thick 
blue line), Galactic proton flux in the escape model (thin 
blue line) together with the resulting photon (red line) and 
neutrino (magenta) fluxes. Upper panel: a p = 2.17 and 
E e sc = 3 x 10 11 eV; Middle panel: a p = 2.10 and E esc = 
3 x 10 11 eV; Lower panel: a p = 2.10 and E esc = 10 14 eV. 
Em ax = 10 19 eV for all three panels. CR proton data from 
KASCADE, KASCADE-Grande [5| and Auger (black error- 
bars) p,Js3| . IGRB and EGB from Fermi-LAT (red error- 
bars) [16| and high-energy neutrino flux from IceCube (ma¬ 
genta errorbars) [171] . 


One can see in the upper panel that for a p = 2.17 
and E esc = 3 x 10 11 eV, the photon flux from the BL 
Lac/FR I populations may explain the entire extragalac- 
tic 7 -ray background. This choice of parameters would 
imply that the main part of the observed TeV y-ray is 
of hadronic origin. The synchrotron peak observed in 
the spectra of BL Lacs at lower energies is caused in this 
picture by electrons which can escape from the central 
engine and radiate most of their energy in the weaker 
magnetic field of the surrounding host galaxy. Note that 
the agreement of the observed and the predicted y-ray 
flux is non-trivial, because the model parameters were 
chosen to fit the UHECR protons, rather than the IGRB 
spectrum. The predicted high-energy neutrino flux from 
these AGNs is about a quarter of the IceCube neutrino 
flux, requiring a Galactic contribution to these neutrinos 
at the level of 75% (dashed magenta line). The diffuse 
y-ray and neutrino fluxes at Earth are due to in-situ pro¬ 
duction and cascade emission during CR propagation. 
We show in the figure both the EGB and the IGRB. One 
should note that as long as the in-situ emission domi¬ 
nates, y-rays (and neutrinos) point back to their sources 
and the y-ray emission is then truely part of the EGB if 
the source is detected, and not part of the IGRB. If ex- 
tragalactic magnetic fields have a negligible impact, the 
distinction would become irrelevant. 

For a p = 2.1 and the same grammage (middle panel), 
the photon and neutrino fluxes are somewhat lower. In¬ 
deed, for a CR harder spectrum and for an extragalactic 
CR flux that satisfies the levels observed at very high en¬ 
ergies, less low energy CRs are present. For a p = 2.1 and 
E esc = 3 x 10 11 eV, the photon flux from the BL Lac/FR 
I populations provides a good fit to the IGRB deduced 
by Fermi-LAT [16} , as can be seen in the middle panel. 

The impact of E esc on the secondary fluxes can be seen 
comparing the middle to the lower panel of Fig. [ 6 j For 
the same slope a p = 2.1, a change in the value of E esc 
affects the resulting diffuse y-ray flux much weaker than 
the high-energy neutrino flux. The latter are produced by 
CRs whose interaction depth is r < 1. In this regime, 
the secondary fluxes scale linearly with the grammage. 
In contrast, the contribution to the diffuse y-ray flux of 
CRs with r 1 is practically not affected by a change 
in E esc . For the parameters chosen in the lower panel, 
one can now explain about ~ 60% of the IceCube flux by 
secondary neutrinos from BL Lacs/FR Is. 

Explaining the extragalactic CR flux within the escape 
model, the IGRB/EGB and a large fraction of IceCube 
neutrinos requires a sufficiently large interaction depth 
at the source (as that chosen in the lower panel). There¬ 
fore our scenario favors CR acceleration close to the black 
hole of the BL Lacs/FR Is, in a region where pp inter¬ 
actions dominate over py interactions, such as the polar 
caps. Let us however note that other sites of CR accel¬ 
eration such as the jets or lobes of BL Lacs/FR Is can 
still explain the IGRB/EGB flux, provided that a p is 
somewhat larger than 2.2. In this case, the contribution 
of these sources to the IceCube neutrino flux would be 
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however reduced. 


C. Neutrinos from BL Lacs and IceCube neutrinos 

A remarkable feature of the secondary spectra is that 
the neutrino flux for our reference parameters can be 
comparable to the flux of astrophysical neutrinos mea¬ 
sured by IceCube na , see especially Fig. [G] (lower panel). 
In Fig. El the dashed magenta lines indicate the Galac¬ 
tic neutrino flux required 1 to match the observed Ice- 
Cube signal: The additional Galactic neutrino contribu¬ 
tion varies between ~ 90% and ~ 30%, depending on 
the value of a p . This is in line with the results from 
Refs. [18], which found a high-energy neutrino flux from 
our Galaxy at the level of ~ 50% of the IceCube flux, 
taking a global Galactic cosmic ray spectrum slope of 2.5 
instead of 2.7 (local flux). 

Also, the slope of our neutrino spectrum is close to the 
one measured by IceCube. The slope of the extragalactic 
neutrino spectrum in the IceCube range is a v ~ a p + 6 ~ 
2.4—2.5 for an injection spectrum with a p — 2.1 —2.2 and 
Kolmogorov turbulence, <5 = 1/3. This is an important 
difference with respect to other models that predict a 
harder spectrum, with slopes ~ 2 . 0 , i.e. similar to the 
slope of the injection spectrum of CRs at the sources. In 
such models, either the parent CRs escape freely from 
the sources or lose all their energy in the sources. In our 
model, CRs diffuse in the source before escaping, which 
results in an additional softening by S of the neutrino 
slope. 

Therefore, our model can explain the entire astrophys¬ 
ical neutrino signal observed by IceCube, both in terms 
of the flux level (Galactic and extragalactic contributions 
of the same order of magnitude) and of the slope. 

The increasing size of the neutrino sample with time 
will allow IceCube to constrain the ratio of the Galac¬ 
tic and extragalactic high-energy neutrino fluxes, study¬ 
ing the anisotropy of their arrival directions. Within 
our model, this flux ratio depends strongly on the slope 
of the extragalactic proton flux and provides therefore 
an important information on the extragalactic UHECR 
sources. 


D. Interactions in the host galaxy and galaxy 
cluster 

We verify in this subsection which impact the host 
galaxy and galaxy cluster have on the diffuse CR proton 
flux that effectively escapes from them, and on the pro¬ 
ductions of secondaries. Our main conclusions are that 


1 We estimate the required Galactic contribution assuming an 

E 2 ' 5 spectrum of the Galactic neutrino flux and normalizing 

the total neutrino flux using the first energy bin of the IceCube 

data at 6 X 10 13 eV. 


only CR protons with energies E < 10 16 eV are confined 
in galaxy clusters and only produce a negligible amount 
of secondary 7 -ray and neutrinos. Therefore, the diffuse 
CR proton flux above ~ 10 16 eV, as well as the diffuse 
7 -ray and neutrino fluxes we computed previously (see 
e.g. Fig. El) are not affected. 

Let us first consider the possibility that CRs residing 
in a kpc scale jet interact with the interstellar medium 
of the AGN host galaxy. The energy loss time t pp for 
the typical ISM density n ~ 1 cm -3 is about 3 x 10 7 yr, 
which is longer than the escape time, see Eq. m , even 
for GeV CRs. Thus CRs in the kpc scale jet escape into 
the interstellar medium of the source host galaxy, rather 
than release their energy inside the jet. 

The density of the intracluster medium spread over 
the Mpc scale of the radio lobes has still lower density 
n ~ 10 -2 — ICC 4 cm" 3 , so that the pp energy loss time 
is comparable to or longer than the age of the Universe. 
CRs residing in the radio lobes then escape into the host 
galaxy cluster of the source, rather than dissipate their 
energy in the lobes. 

Therefore, CRs produced in the AGN jet or in the radio 
lobes escape in the host galaxy and galaxy cluster. For a 
magnetic field strength of B ~ 1 pG, which is typical for 
galaxy clusters within a Mpc region, the escape time of 
the very and ultra-high energy CRs contributing to the 
extragalactic proton flux at Earth (above ~ 10 16 eV) is 
small compared to the age of the Universe. 

Lower energy cosmic rays produced by an UHECR 
source which operates only a limited time (~ 10 ® yr in 
the case of radio-loud AGNs) are still found in the cluster 
long since the UHECR source has ceased to exist. Part 
of their energy will be released while residing in the host 
galaxy and galaxy cluster. Calculating E esc for the case 
of the host galaxy cluster, one finds that relativistic par¬ 
ticles do not escape for n ~ 10~ 4 cm _3 . Host galaxies 
and galaxy clusters are therefore expected to give only 
a subdominant contribution to the diffuse neutrino flux. 
As discussed above, the diffuse 7 -ray flux is less sensi¬ 
tive to the value of E esc and, pp interactions in the host 
galaxy may contribute to the IGRB depending on the 
slope a p . 


V. CONCLUSIONS 

We have found that star-forming galaxies (normal spi¬ 
ral galaxies, starburst galaxies and SF-AGNs) give only 
a sub-dominant contribution to the high-energy CR flux 
(E < 10 18 eV). Both their overall luminosity, their red- 
shift evolution which peaks at z ~ 2 — 3, and their 
maximal energy disfavor this source class as the main 
source of extragalactic CRs up to the ankle. The result¬ 
ing secondary fluxes are more uncertain, since they de¬ 
pend on the not well determined grammage CRs cross in 
their host galaxies and galaxy clusters. Even keeping the 
grammage as a free parameter, it is not possible to ex¬ 
plain at the same time a large contribution to the IGRB 
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and to the IceCube neutrinos or to the extragalactic part 
of the IceCube neutrino signal. 

In contrast, we have shown that the BL Lac/FR I pop¬ 
ulation as a source for extragalactic CRs can explain in 
a unified way both the observations of primary and sec¬ 
ondary fluxes. The main reason for this difference with 
star-forming galaxies or other sources is that the num¬ 
ber density of BL Lac and FR I galaxies peaks at low 
redshift, z < 1 . 

More precisely, we found that the extragalactic CR 
proton flux can be explained for any acceleration site 
(close to the black hole, in the jets or in the lobes) of 
the BL Lacs/FR I galaxies. However, only acceleration 
close to the black hole (especially at the polar caps, see 
Section III Bh satisfies the required conditions to produce 
secondary 7 -ray and neutrino fluxes that can explain the 
extragalactic IceCube neutrino flux and of the IGRB. 
For a spectral slope of CR protons close to a p = 2.2, as 
suggested by shock acceleration and the escape model, 
we find that such UHECR sources provide the dominant 
fraction of both the isotropic 7 -ray background and of 
the extragalactic part of the astrophysical neutrino signal 
observed by IceCube. 

We showed that the difference in the slopes of the pro¬ 


ton and the neutrino fluxes can be explained by the dif¬ 
fusion of primary protons in the turbulent magnetic field 
of CR sources. In the case of Kolmogorov turbulence, the 
power-law of the secondary neutrino spectrum is changed 
by 1/3, what explains the relatively soft neutrino spec¬ 
trum with a„ ~ 2.5 observed by IceCube using a proton 
spectrum with a p — 2.1 — 2.2. This mechanism is uni¬ 
versal and does not depend on the type of sources. 
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